ApJ Letters, in press 

Preprint typeset using L£-T^X style cmulateapj v. 5/2/11 



THE PROTOCLUSTER G18.67+0.03: A TEST CASE FOR CLASS I CH 3 OH MASERS AS EVOLUTIONARY 
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ABSTRACT 



R. Indebetouw • , C.J. Chandler" 



We present high angular resolution Submillimeter Array (SMA) and Karl G. Jansky Very Large 
Array (VLA) observations of the massive protocluster G18.67+0.03. Previously targeted in maser 
surveys of GLIMPSE Extended Green Objects (EGOs), this cluster contains three Class I CH3OH 
maser sources, providing a unique opportunity to test the proposed role of Class I masers as evolu- 
tionary indicators for massive star formation. The millimeter observations reveal bipolar molecular 
outflows, traced by 13 CO(2-l) emission, associated with all three Class I maser sources. Two of these 
sources (including the EGO) are also associated with 6.7 GHz Class II CH3OH masers; the Class II 
masers are coincident with millimeter continuum cores that exhibit hot core line emission and drive 
active outflows, as indicated by the detection of SiO(5-4). In these cases, the Class I masers are 
coincident with outflow lobes, and appear as clear cases of excitation by active outflows. In contrast, 
the third Class I source is associated with an ultracompact HII region, and not with Class II masers. 
The lack of SiO emission suggests the 13 CO outflow is a relic, consistent with its longer dynamical 
timescale. Our data show that massive young stellar objects associated only with Class I masers are 
not necessarily young, and provide the first unambiguous evidence that Class I masers may be excited 
by both young (hot core) and older (UC HII) MYSOs within the same protocluster. 

Subject headings: ISM: individual objects (GI8.67+0.03) — ISM: jets and outflows — ISM: molecules 
— masers — stars: formation — techniques: interferometric 



1. INTRODUCTION 

The lack of a detailed, observationally based evolution- 
ary sequence for massive young stellar objects (MYSOs) 
limits our understanding of the early stages of high mass 
(Mzams >8 M Q ) star formation. Massive stars form in 
(proto)clusters, with younger sources (such as hot cores) 
often found in c lose proximity to ultracompact (UC) 
HII r egions (e.g. iHunter et al.l 120061 : iCvganowski et all 
120071 and references therein). Because massive proto- 
clusters are generally distant (D>I kpc) and deeply em- 
bedded, most efforts to develop evolutionary sequences 
have focused on cm-wavelength maser transitions, which 
are amenable to high resolution study and unaffected 
by extinction. Four types of masers are ubiquitous: 
Class I and II CH 3 OH, H 2 0, and OH. Class I CH 3 OH 
and H2O masers are collisionally pumped in shocked 
gas, while Class II CH3OH and OH masers are radia- 
tively pumped by infrared emission from warm dust (e.g. 
Elitzur et ail 119891: iCragg et al.l 120021: iVoronkov etaLl 
20061 ). CHqOH masers are key to proposed evolutionary 
sequences of masers in MYSOs, w hich posit that Class 
I CH3OH masers a ppear first (e.g. lEllingsen et al.ll2007t 
iBreen et al.ll2010al ). based in part on their observed as- 
sociation with outflows (e.g. iPlambeck fc Mentenl 119901 : 
iKurtz et al J 12004 : ICvganowski et all 120091) . Class I-onlv 
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sources have thus generally been interpreted as tracing, 
the e arliest stages of massive star formation (|Ellingsenl 
2006). Recent work, however, suggests that Class I 
CH3OH masers may also be excited in shocks driven 
by expanding HII regions (jVoronkov et al.ll2010l ). rais- 
ing the possibility of multiple distinct epo chs of Class 

I maser activity during M YSO evolution (jChen et al.l 
l20Tlt IVoronkov etai] [2012]) . The prevalence of Class I- 
only sources is currently unknown, primarily due to the 
lack of substantial untargeted searches in Class I maser 
transitions (see also IVoronkov et al.ll2010l ). 

Our studies of a new sample of MYSO outflow 
candidates (Extended Green Objects (EGOs), selected 
based on extended 4.5 / xm em ission in GLIMPSE im- 
ages; ICvganowski et all I2008D have ide ntified a mas- 
sive protocluster (M c i ump ~3510 M ; iSchuller et al.1 
120091 ) that provides a unique opportunity to un- 
derstand the evolution of Cla ss I CH3OH masers: 

G18.67+0.03. ICvganowski eTaTl (|200flh detected 44 GHz 
Class I CH3OH masers towards three mid-infrared (MIR) 
sources in this region: two (including the EGO) are 
also associated with 6.7 GHz Class II CH3OH masers 
(Fig. [T]). The third-the only source in the region with 
a cm- wavelength continuum counterpart in deep Ver y 
Large Array (VLA) images (ICvg anows ki et aT1l2011b[)- 
has no Class II masers. iGreen fc McClure-Gr iffiths 
(|2011l ) assign both 6.7 GHz masers in G18.67+0.03 the 
far distance in their HI self-absorption study of Class 

II masers. Using the r evised kinematic distance pre- 
scription of iReid et al.1 ( 20091 ) and thermal NH 3 data, 
ICvganowski et al. I ( 2012 ) find a similar far distance of 
10.8 kpc, which we adopt here. 

In this Letter, we present a high-resolution Submillime- 
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Figure 1. Three color Spitzer image (RGB: 8.0,4.5,3.6 jum) of G18.67+0.03; positions of CH 3 OH masers from ICvganowski et al.l II2009D are 
marked with O (Class II) and + (Class I). In all panels, the synthesized bcam(s) are shown at lower left, (a) SMA pointings (dashed yellow circles) 
and contours of BGPS 1.1 mm emission (white) and VLA 3.6 cm emission (blue; Cyganowski ct al. 2011b) arc overlaid. Contour levels: BGPS: 
(10,20) x<7, a =0.03 Jy beam -1 ; VLA: (4,40,120) X a, a =3.06E-5 Jy beam -1 , (b) Contours of SMA 1.3 mm continuum emission (gray) and MIPS 24 
(im emission (yellow). Levels: SMA: (3,5,10,20,40,60,120) xcr, a =0.00145 Jy beam -1 . MIPS: 600,1200,1800 MJy sr -1 . (c-c): Zoom views of maser 
sources, with SMA contours (black) as in (b). In (d), VLA 1.2 cm continuum contours (blue) are also overlaid (levels (5,10,25,50) xcr, cr =7.4E-5 Jy 
beam -1 ). 

ter Array (SMA j3 study of the mm-wavelength dust con- 
tinuum and molecular line emission toward G18.67+0.03, 
to constrain the evolutionary states of the maser sources 
and the presence/absence of active outflows. The SMA 
is well-suited for deciphering MYSO relative evolu- 
tionary states (jBeuther et all 120091: I Zhang et all 120091: 



7 The Submillimeter Array is a joint project between the Smith- 
sonian Astrophysical Observatory and the Academia Sinica Insti- 
tute of Astronomy and Astrophysics and is funded by the Smith- 
sonian Institution and the Academia Sinica. 



iWang et alll2012T ) . To provide kinematic information for 
the ionized gas, we include complementary 1.3 cm Karl 
G. Jansky Very Large Array (VLA) hydrogen recombi- 
nation line and cm continuum observations. 

2. OBSERVATIONS 

SMA observations at 1.3mm were obtained on 18 June 
2010 with 7 antennas in compact configuration, in good 
weather (jii^ghz ~0.07). A three-pointing mosaic was 
used to cover the extent of the (sub)mm clump (Fig.[T]). 
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The primary beam (FWHP) at 1.3 mm is ~55". The 
largest angular scale recoverable from these data is ~20". 
We observed ^216.8-220.8 GHz in the lower sideband 
(LSB) and -228.8-232.8 GHz in the upper sideband 
(USB), each divided into 2x2 GHz IFs, with a uni- 
form channel width of 0.8125 MHz. The calibrators were 
J1733-130and J1911-201 (complex gain) , 3C454.3 (band- 
pass), and Neptune (absolute flux density). The LSB 
was flux-calibrated using a model of Neptune's bright- 
ness distribution and the MIRIAD task smaflux. This 
direct approach cannot be used for the USB because of 
strong 12 CO in Neptune's atmosphere. Instead, the USB 
flux density calibration was bootstrapped from the LSB 
assuming a spectral index a = —0.85 for 3C454.3 (based 
on SMA monitoring) . Comparison of derived quasar flux 
densities with SMA monitoring suggests the absolute flux 
density calibration is accurate to better than ~15%. 

Initial calibration was performed in MIRIAD; the data 
were then exported to CASA for further processing. 
Each IF was processed independently, with the contin- 
uum estimated using line-free channels in the Mw-plane, 
and separated from the line emission. The continuum 
was self-calibrated and the resulting solutions applied 
to the line data. The data were imaged using Briggs 
weighting and a robust parameter of 0.5. The final com- 
bined continuum image has a synthesized beamsize of 
3"16x2"95 (P.A.=21°) and a lcr rms of 1.5 mJy beam" 1 . 
The line data were resampled to Av=1.2 km s -1 , then 
Hanning-smoothed; the typical rms per channel is ~40 
mJy beam -1 . All measurements were made from images 
corrected for the primary beam response. 

The NRAC0 VLA observations of G18.67+0.03 were 
obtained in the C configuration on 7 January 2011 as part 
of a 1.3 cm line and continuum su rvey of massiye pro - 
tostellar objects (program AB1346; iBrogan et alj |2011). 
The calibrators were J1832-1035 (complex gain), J1924- 
2914 (bandpass), and 3C286 (absolute flux density). A 
single pointing was observed (a =18 h 24 m 53 s .700, 8 = 
-12°39'20'.'0); the primary beam (FWHP) is -2'. The 
data were calibrated and imaged in CASA. The full VLA 
dataset will be described in a future publication; here we 
consider only the 24.50878 GHz H64a, 25.68788 GHz 
H63a, and 24.64 GHz (1.2 cm) continuum data. The 
self-calibrated continuum image has a synthesized beam- 
size of T/46xT/08 (P.A.= -174°) and a la rms of 74 
/iJy beam -1 . To improve S/N, we smoothed the H63a 
and H64a data (to Av channcl =4.8 km s -1 , 6 beam =2"), 
then averaged to obtain a final radio recombination line 
(RRL) image cube with a ^0.29 mJy beam -1 . All mea- 
surements were made from images corrected for the pri- 
mary beam response. 

3. RESULTS 

3.1. Millimeter Continuum and Line Emission 

As shown in Figure [TJ we resolve four compact 1.3 
mm continuum sources, arranged in a line of ~50" (2.6 
pc) extent along the clump. All have 24 /xm counter- 
parts. For ease of reference, we designate them (from 
E to W): EGO, CM1-UCHII, Non-EGO maser source 

8 The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under agreement by the 
Associated Universities, Inc. 



(NEMS), and mid-infrared source (MIRS). CM1-UCHII 
refers to the mm counterpart of the cm-wavelength 
continuum source design ated F G18.67+0.03-CM1 by 
ICvganowski et all (|2011b|) . Observed millimeter contin- 
uum properties are summarized in Table [TJ The SMA 
image recovers 19±6% of the 1.1 mm BGPS flux den- 
sity (3.8t° 8 7 Jy; iRosolowskv et al.|[2010l corrected as per 
lAguirre et al.l [20Tll) . The westernmost source (MIRS) 
fal ls outside the area searche d for 44 GHz CH3OH masers 
by ICvganowski et al.l ([2009;), so we do not discuss it fur- 
ther in this Letter. 

Only a few species (including CO isotopes, SO, DCN, 
and low-excitation lines of H2CO and CH3OH) are de- 
tected towards all three maser sources. The LSR veloc- 
ities of the sources, determined from compact molecular 
line emission, are within 3 km s -1 of each other (Fig.[3Ji, 
Table [JJ . Figure [3Ja,-c presents a comparison of the spec- 
tra at the EGO, CM1-UCHII, and NEMS 1.3 mm con- 
tinuum peaks across 2 GHz of the SMA band. The EGO 
and NEMS both exhibit copious molecular line emission 
from classic hot core tracers, such as CH3CN, OCS, and 
high-excitation CH 3 OH. In contrast, CM1-UCHII is line- 
poor: none of these tracers are detected. 

All three Class I CH3OH maser sources are associated 
with both redshifted and blueshifted 13 CO(2-l) emission, 
with velocity extents of 23, 16, and 30 km s -1 for the 
EGO, CM1-UCHII, and NEMS, respectively (Fig. [3b). 
In all cases, there is a clear spatial offset between the 
red and blue lobes, which are centered on a compact 
source. These characteristics indicate that the high- 
velocity 13 CO(2-l) emission traces bipolar molecular out- 
flows. The 44 GHz masers associated with the EGO 
and NEMS are coincident with outflow lobes. The CM1- 
UCHII 44 GHz masers are located, in projection, near the 
edge of its redshifted lobe. In addition to the N-S out- 
flow, there is evidence for an E-W gradient in low- velocity 
molecular gas near the UCHII region (v— v sys t e mic J$L5 
km s -1 ; Fig. [5b). This E-W velocity gradient could in- 
dicate a second, slower outflow, and/or large-scale rota- 
tion or infall. Importantly, the outflows associated with 
the EGO and NEMS are detected in SiO(5-4) emission, 
while the outflow(s) associated with CMI-UOHII is not 
(Fig. HJi). 

3.2. Temperatures and Dust Masses 

We estimate gas temperatures for the EGO and 
NEMS by fitting the J =12-ll CH 3 CN ladder u sing 
the method described in ICvganowski etldl (|2011aD . In 
brief, the model accounts for optical depth effects and 
emission from the CHg 3 CN isotope: the temperature, 
source size(diameter), and CH3CN column density and 
linewidth are free parameters. Figure[5Ji-e show the best- 
fit models, overlaid on the observed spectra. The best-fit 
temperatures, source sizes, and CH3CN column densities 
are quite similar for the two sources: T=175K and 185K, 
log[N CH3 CN(cm -2 )]=16.66 and 16.58, and d source =0'/45 
(4900 AU) and 0'.'39 (4200 AU) for the EGO and NEMS, 
respectively. 

Gas masses and (beam-averaged) column and num- 
ber densities estimated from the 1.3 mm thermal dust 
emissio n, using the methodology of Cvganowsk i et all 
(|2011al) . are presented in Table [TJ The mass estimates for 
the EGO and NEMS assuming T dust =T C H 3 CN (~30 and 
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Table 1 

Properties of Millimeter Continuum Sources 





Observed Properties 








Derived Properties 




Sourcc a 


J2000 Coordinates 15 


Peak Intensity Intcg. 


Flux Density T du3t 


Tduat 


M ga3 Njj 2 d n H2 d 


VLSR° 


Name 


Q( h m B) g(0 , „ } 


(mjy beam -1 ) 


(mjy) 


(K) 




(M ) (10 23 cm" 2 ) (10 5 cm" 


3 ) km s" 1 


EGO 


18 24 53.77 -12 39 20.8 


96 


153(36) 


175 


0.003 


30 0.7 2.8 


79 


CMl-UCHII f 


18 24 52.59 -12 39 20.0 


36 


126(41) 


50 


0.008 


100 2.1 8.6 


80 


NEMS 


18 24 51.09 -12 39 22.0 


187 


347(44) 


185 


0.006 


70 1.5 6.1 


82 


MIRS 


18 24 50.27 -12 39 22.0 


44 


92(32) 








81 



Total 718 



a As designated in EI3.1I 

Coordinates of 1.3 mm continuum peak. The number of significant figures reflects a one-pixel uncertainty. 
c Integrated flux density within 3(7 contour, measured using the CASAVIEWER program. Following the approach developed 
for the CORNISH survey iPurcell et al.ll2008D . we estimate uncertainties as ^/./V sr cr 2 ky (l 4- ) , where N Br is the number of 

pixels in the source aperture, N s k y is the number of pixels in an off-source "sky" aperture, and cr s k y is the standard deviation 
over the "sky" aperture. The quoted uncertainties are averages for 3 choices of "sky" apertures. 
d Beam-averaged 

c Determined from compact molecular line emission, i|3. II 

* Physical properties calculated using the 1.3 mm flux density less the estimated free-free contribution of 5 mjy. 
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Figure 2. SMA spectra towards the (a) EGO, (b) CM1-UCHII, and (c) NEMS 1.3 mm continuum peaks, showing 
(d-e) SMA CH3CN spectra (gray) overlaid with the best- fit model (dashed black line). 



2 GHz of the LSB. 



70 M Q , respectively) are likely lower limits; the CH3CN- 
emitting regions are unresolved at the scale of our SMA 
observations (9 syn . ~33,000 AU), while the 1.3 mm con- 
tinuum emission appears somewhat more extended. For 
comparison, the gas masses of the CHaCN-emitting re- 
gions, estimated from the best-fit source sizes and column 
densities, are ~40 and 30 M for the EGO and NEMS, 
for CH 3 CN/H 2 =10~ 8 . Since the CH 3 CN/H 2 abundance 
in hot cores is uncertain, and values an order of magni- 
tude h igher /lower are plausible (see ICveanowski et al] 
1201 laf) . this is reasonable agreement. 

Constraining the gas temperature of CM1-UCHII is 
more difficult because of the paucity of associated line 



emission: no CH3CN is detected, and only one CH3OH 
line. Most lines detected have E upper <50 K; the highest- 
excitation transitions have E uppcr ^68 K (H2CO). In ad- 



dition, while the ionized gas is unresolved f c]3.3[) . the mm- 
wavelength continuum emission is extended E-W, with 
scale >0.5 pc. For the estimates in Table [TJ we adopt a 
temperature of 50 K, and subtract the free- free contribu- 
tion to the 1.3 mm flux density (extrapolated from the 
1.2 cm integrated flux density assuming optically thin 
emission, S y oc v~ 
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Figure 3. SMA molecular gas kinematics, (a) DCN(3-2) first moment map. (b) Zoom view of CM1-UCHII H 2 CO(3o,3-2o,2) velocity field 
(moment 1, colorscale), overlaid with contours of VLA 1.2 cm continuum emission (black) and blue/red-shifted 13 CO(2-l) emission, (c) SMA 1.3 
mm continuum (colorscale) overlaid with contours of VLA 1.2 cm continuum (black) and composite maps of blue/red-shifted 13 CO(2-l) emission. 

CO is integrated over the indicated velocity ranges, optimized for each source. Contour levels (Jy beam" - 1 km s" 1 ): EGO: blue: 1,2,3,4,5 : red: 
1,2,3,4. CM1-UCHII: blue: 2,3; red: 2,3,5,7. NEMS: blue: 2,4,6; red: 2,4,6,8. MIRS: blue: 1.5,2,3,4; red: 1.5,2,3,4. (d) SiO(5-4) integrated intensity 
map (moment 0, colorscale), o verlaid with VLA 1.2 cm continuum contours (white). VLA contour levels as in Figure^ Positions of CH3OH masers 
from ICvganowski ct al. ||2009D are marked with O (Class II) and + (Class I). Synthesized bcam(s) arc shown at lower left in each panel. 

3.3. Ionized Gas Properties 

We detect the H63a and H64a recombination line s 
towards G18.67+0.03-CM1 (ICvganowski et all l20TTbh . 



The fitted RRL velocity (77.3±0.9 km s" 1 ) is consis- 
tent with the molecular-gas velocity. From the line-to- 
continuum ratio of 0.48, the fitted FWHM linewidth of 
27.0 km s _1 , and the fitted source size of 0'.'63 (from 
the continuum image), we derive an electron tempera- 
ture of ^5100 K and de nsity of ~1.9xl0 4 c m~ 3 (follow- 
ing the method of e.g. iGarav et al.l 119861 : iSewilo et al.l 



[20T1 . If bulk motions (e.g. expansion) dominate the 
linewidth, this could explain the low derived T e ; however, 
the RRL emission is too weak to investigate this possibil- 
ity. The cm-wavelength spectral index is consistent with 
optically thin free-free emission, and t he ionizing photon 
flux (~ lxlQ 47 s _1 , calculated as in ICvganowski et al.l 
1201 lbf) corre sponds to a single ionizing star of spectral 
type B0.5V (ISmith et al.H2002D . 
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4. DISCUSSION 

4.1. Evolutionary States of Class I Maser Driving 
Sources 

All three maser sources in G18.67+0.03 are asso- 
ciated with massive (>30 M Q ), dense (n Hs , >10 5 
cm -3 ) millimeter continuum cores, and with 13 CO out- 
flows. CO (a nd HCO + ) may trace relic outflows from 
MYSOs (e.g. IKlaassen et all [20061 : IKlaassen fc Wilsonl 
120071: IHunter et al.l I2008D . SiO emission, however, 
specifically indicates recent shocks and active out- 
flows, since the gas-phase SiO abundance is enhanced 
by shocks and remains so f or only ~10 4 years (e.g. 
IPineau des Forets et al.l (1997?) . The EGO and NEMS, 
which have both Class I and II CH3OH masers, thus 
share two other key characteristics: (1) their outflows 
are active, as evidenced by SiO(5-4) emission; and (2) 
the driving sources exhibit hot core line emission (e.g. 
CH3CN, CH3OH, OCS). In summary, the Class I masers 
in these sources appear as clear cases of excitation by 
outflows. 

In contrast, the Class I-only maser source lacks indica- 
tors of youth: neither SiO nor hot core line emission is de- 
tected in our observations. The 44 GHz masers lie near a 
cm continuum source, detected in RRLs; its derived ion- 
ized gas properties are consistent with a UCffll region. 
Two mechanisms have been posited that could explain 
Class I maser emission associated with such a compara- 
tively evolved source. IVoronkov et al.l (|2006f ) suggested 
that Class I masers might be long-lived, persisting af- 
ter the exciting shoc k had dissipated. More recently, 
IVoronkov et al.l ( |2010l ) noted the proximity of cm contin- 
uum emission to three of the four known examples of 9.9 
GHz Class I CH3OH masers. They suggested that shocks 
driven into the surrounding cloud by expanding HII re- 
gions were responsible for these masers, and that this 
mechanism should also apply to other Class I transitions. 
These two scenarios have somewhat different implica- 
tions for maser evolutionary sequences: the first implies 
that Class I masers outlast the Class II phase; the second, 
that Class I masers may appear more th an once, excited 
by yo ung or evolved MYSOs (see also IVoronkov et al.l 
l2012h . 

In G18.67+0.03, the lack of associated SiO emission 
suggests that the CM1-UCHII 13 CO outflow is a relic. 
It also has a smaller velocity extent than the EGO and 
NEMS 13 CO outflows, despite being more spatially ex- 
tended: additional evidence that the CM1-UCHII out- 
flow is older (t dyn ~ 19000, 26000, and 66000 for the 
NEMS, EGO, and CM1-UCHII respectively). In this 
picture, the CM1-UCHII 13 CO outflow would have been 
driven by the now-ionizing source, prior to the creation 
of the UC HII region (analogou s to the E-W outflow in 
G5.89-0.39: DHunter et al.ll200l . The Class I masers are 
near the edge of the redshifted lobe, consistent with an 
association with the relic(?) flow, and so supporting the 
scenario in which long-lasting Class I maser activity ex- 
tends beyond the Class II maser lifetime. This interpre- 
tation, however, is not conclusive; the Class I masers are 
also near the UCHII region, and we cannot rule out their 
being excited by HII region-driven shocks. In either case, 
these Class I masers are clearly associated with a more 
evolved MYSO. 




60 70 80 90 100 110 120 

Velocity (km s _1 ) 




60 70 80 90 100 110 120 

Velocity (km s _1 ) 



Figure 4. SMA CH3OH spectra at the position of the strongest 44 
GHz maser associated with (a) NEMS and (b) CM1-UCHII. Solid 
line: 8_i 8 -7o 7 (229.759 GHz); dashed line: 3_ 2 2-4_i 4 (230.027 
GHz); dotted line: 8 ,8-7i, 6 (220.078 GHz). 

4.2. Class I Maser Excitation as a Key? 

Probable maser emission in the Class I 229.759 GHz 
CH30H(8_i i 8-7o,7)E transition is a conspicuous feature 
of recent SMA observat ions of MYSO outflows (e.g. 
I Cyganowski et afl l2011al and references therein). In 
these outflows, 229 GHz emission is generally colocated 
(spatially and spectrally) with lower frequency Class I 
masers, and the brightest featur es in all transitions co- 
incide (see also iFish et "all 1201 11 ) . In G18.67+0.03, we 
detect strong 229.759 GHz emission coincident with the 
brightest NEMS 44 GHz masers (Fig. SJ-those associ- 
ated with the blueshifted, SiO-rich western outflow lobe. 
Like many previous 229 GHz studies, our beam is too 
large to establish masing based on brightness tempera- 
ture. iSlvsh et al.l (|2002T ) proposed the 229.759/230.027 
ratio as a diagnostic, with values >3 indicating nonther- 
mal 229 GHz emission: the ratio for the NEMS position 
shown in Figure|4]is 5.5. Thi s NEMS 44 GHz maser is the 
second-brightest detected bv lCvganowski et al.1 (|2009f ) in 
G18.67+0.03; the brightest is the western of the two as- 
sociated with CM1-UCHII (Fig. [TJi; I pC ak,44GHz ~2 Jy 
beam -1 ). No 229.759 GHz emission is detected at this 
position (Fig. 0]). The 44 GHz masers associated with 
the EGO are weak (< 0.5 Jy beam -1 ), and the mm line 
ratios are ambiguous. 

These results suggest that excitation conditions may 
differ in older and younger Class I maser sources in a 
manner that produces observable differences in the rela- 
tive strengths of various Class I maser transitions. De- 
tailed comparisons of the 44 and 229 GHz maser spots 
are, however, limited by the relatively coarse spatial and 
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spectral resolution of the current SMA data. In gen- 
eral, additional modeling and observations are needed to 
understand what information about physical conditions 
and evolutionary state can be extracted from the pres- 
ence/absence of different Class I CH3OH masers, particu- 
lar ly for mm-wavel ength transitions (see also discussion 
in iFish et al1l20lTI ). For G18.67+0.03, higher angular 
resolution observations are needed to confirm the maser 
nature of the 229 GHz emission, and to spatially resolve 
the potential mm masers from the hot cores in the EGO 
and NEMS; high r esolution 36 GHz obs ervations may 
also be helpful (e.g. iVoronkov et aDl2012| ). 

5. SUMMARY 

Our high resolution SMA and VLA observations pro- 
vide the first unambiguous evidence of Class I CH3OH 
masers being excited by both young and more evolved 
MYSOs within the same protocluster. The two hot cores, 
also associated with Class II CH3OH masers, drive ac- 
tive outflows traced by SiO(5-4) emission and Class I 
masers. In contrast, the UCHII region is associated 
only with Class I CH3OH masers and with an older 
(possibly relic) 13 CO outflow; the UCHII shows nei- 
ther SiO nor hot core line emission. These results fur- 
ther demonstrate the limitations of curr ent evolutionary 
sequences for maser em ission (see also IVoronkov et al.l 
2010; iChen et~aT|[20TTI) . In particular, our data show 
that MYSOs associated only with Class I CH3OH masers 
(no Class II) are not necessarily young, in contrast to 
published sequences. This work highlights the impor- 
tance of high angular resolution multiwavelength obser- 
vations for constraining MYSO evolutionary states, and 
disentangling an observation-based MYSO evolutionary 
sequence. 

This research made use of NASA's Astrophysics 
Data System Bibliographic Services, APLpy, an 
open-source plotting package for Python hosted at 
http://aplpy.github.com, and the myXCLASS program 
(http://www.astro.uni-koeln.de/projects/schilke/XCLASS), 
which access es the CDMS (|http://www.cdms.de| and 
JPL (http://spec.jpl.nasa.gov) molecular databases. 
C.J. Cyganowski is supported by an NSF AAPF under 
award AST-1003134. 
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